Bose-Einstein Condensation of S = 1 Ni spin degrees of freedom in NiCl2-4SC(NH 2 ) 2 
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It has recently been suggested that the organic compound NiCl2-4SC(NH2)2 (DTN) exhibits 
Bose-Einstein Condensation (BEC) of the Ni spin degrees of freedom for fields applied along the 
tetragonal c-axis. The Ni spins exhibit 3D XY-type antiferromagnetic order above a field-induced 
quantum critical point at H c i ~ 2 T. The Ni spin fluid can be characterized as a system of effective 
bosons with a hard-core repulsive interaction in which the antiferromagnetic state corresponds to a 
Bose-Einstein condensate (BEC) of the phase coherent S = 1 Ni spin system. We have investigated 
the the high-field phase diagram and the occurrence of BEC in DTN by means of specific heat and 
magnetocaloric effect measurements to dilution refrigerator temperatures. Our results indicate that 
a key prediction of BEC is satisfied; the magnetic field-temperature quantum phase transition line 
H C (T) — H c i oc T a approaches a power-law at low temperatures, with an exponent a — 1.47 ± 0.06 
at the quantum critical point, consistent with the BEC theory prediction of a = 1.5. 

PACS numbers: 75.40.-s,65.40.Ba, 
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In the past few years, a class of quantum spin sys- 
tems that exhibit a Bose-Einstein condensation in ap- 
plied magnetic fields has been receiving an increasing 
amount of attention. The compounds studied previously 
consist of chains of Ni atoms planes of Cu dimcrs, 
(BaCuSi0 6 )i£ or 3-D coupled spin ladders (TlCuCl 3 and 
KCuCl3)pii£ In these compounds, the spin degrees of free- 
dom are generally provided by either weakly interacting 
S — 1 Ni chains exhibiting a Haldane gap, or by a lat- 
tice of dimerized 5 = 1/2 Cu 2+ ions. The spin singlet 
ground state of these systems is separated from the low- 
est excited state (triplet) by a finite energy gap at zero 
magnetic field. In the presence of a magnetic field, the 
Zeeman term reduces the energy of the S = — 1 triplet 
state, until it reaches the energy of the non-magnetic 
ground state at a field H&. Canted XY antifcrromag- 
netism is then observed between H c \ and an upper criti- 
cal field H C 2. This XY antiferromagnetic phase of S = 1 
spins has a U(l) symmetry and can thus be interpreted as 
a BEC. The U(l) symmetry requires a spin environment 
that is axially symmetric with respect to the applied field. 
Observation of a gapless Goldstone mode (magnons) in 
the ordered phase could in principle provide evidence of 
a broken continuous U(l) symmetry. Inelastic neutron 
scattering measurements of the compound TICUCI3 were 
proposed to be consistent with this interpretation^ How- 
ever, ESR measurements subsequently revealed signifi- 
cant magnetic anisotropy in this compound inconsistent 
with a U(l) symmetry of the spins£ In general, axial 
symmetry is not expected for weakly coupled spin ladder 
and spin chain systems. 

The compound NiCl 2 -4SC(NH 2 )2 (DTN)2 is a new 
candidate for Bose-Einstein condensation of spins, and 
has several features that make it unique among these 
BEC of spin systems. In a similar manner to BaCuSi206, 
the tetragonal crystal symmetry satisfies the U(l) sym- 
metry requirement for BEC with the additional feature 



that the symmetry can be tuned by rotating the applied 
magnetic field. It has been predicted^ that an XY mag- 
net should occur for fields along the c-axis, and an Ising 
magnet for angles up to 40 degrees away from the c- 
axis. The spin configuration is also different from previ- 
ous BEC compounds, consisting only of a spin-one state 
with no spin-singlet state. The tetragonal spin lattice 
provided by the Ni ions is described by the Hamiltonian: 



(1) 



where v = {a,b,c} and J a = Jt, due to the tetrago- 
nal symmetry of the crystal structure. The single ion 
anisotropy D ~ 10 Kl2iil splits the the Ni S = 1 spin 
state into the S z = ground state and the S z = ±1 ex- 
cited states. Thus, for an applied magnetic field along 
the c-axis, the level crossing occurs between two triplet 
states rather than a triplet and a singlet as is the case 
for all other BEC systems studied to date. Magnetiza- 
tion measurements^ have revealed AFM order between 
H c i ~ 2 T and H C 2 ~ 12 T with a maximum Neel tem- 
perature of 1.2 K. For a magnetic field perpendicular to 
the c axis, no transition is observed up to H = 15 T 
since such a field mixes the S z = state with a linear 
combination of the S z = ±1 states, producing an effec- 
tive repulsion between energy levels that precludes any 
field-induced quantum phase transition. 

In this work, we investigate the temperature-magnetic 
field phase diagram for H | | c via thermodynamic mea- 
surements. A key prediction of the BEC theory is a 
power-law temperature dependence of the phase transi- 
tion line H C (T) — H c \ oc T™ where a = 1.5^ Previous 
studies on DTN and other candidate BEC spin systems 
have found a wide variety of values of a, generally higher 
than the predicted value of 1.5. This is due to the fact 
that the fits are performed at relatively high tempera- 
tures, whereas the power-law universal behavior is only 
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FIG. 1: Specific heat C versus temperature T of N1CI2- 
4SC(NH 2 ) 2 for H \\ c in applied magnetic fields of 0, 2, 3.5, 
and 10 T. 
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FIG. 2: Magnetocaloric effect data determined by monitoring 
T while sweeping B up and down. Solid arrows indicate the 
direction of B and dashed arrows indicate the inflection point 
in T versus B. 



valid for T c — > 0£ Furthermore, the value of a is very 
dependent on the value of H c i chosen for the fit, which is 
difficult to extrapolate from higher temperature data. A 
previous work found a = 2.6 from fits to magnetization 
data on DTN with most of the data lying between 0.5 and 
1.5 In this presetn study, we obtain a from specific 
heat and magnetocaloric effect measurements to dilution 
refrigerator temperatures. We determine the value of a 
using the extrapolation method described by Sebastian 
et al& The validity of this extrapolation has been verified 
using Monte Carlo calculations of H C (T) and the method 
has been successfully used to obtain the exponent for the 
compound BaCuSi206p£ 

Specific heat and the magnetocaloric effect (MCE) 
were measured in a dilution refrigerator system with 



an 18 T magnet at the National High Magnetic Field 
Laboratory in Los Alamos, NM. All measurements were 
conducted on single crystals with the external field ori- 
ented along the tetragonal c-axis. Specific heat (100 mK 
< T < 1.5 K) was determined by the quasi-adiabatic 
heat pulse relaxation method, and the magnetocaloric 
effect was measured by sweeping the field up an down at 
0.1 T/min while monitoring the temperature between 80 
mK and 1.5 K. 

Distinct thermodynamic transitions can be observed in 
the specific heat and magnetocaloric effect data at mag- 
netic fields between 2.1 and 12.7 T. Representative data is 
shown in Figs ^ and [5] and the phase diagram determined 
from these data is shown in Fig. [31 The specific heat 
data shown in Fig. ^exhibits sharp peaks for H = 3.5 T 
and 10 T. These peaks resemble the A-like transition of 
the specific heat in superfluid helium, an archetypal BEC 
system. Detailed theoretical predictions for this system 
are necessary to make a precise comparison. An equal en- 
tropy construction was used to determine the midpoint of 
the transition. For transitions occurring at a given tem- 
perature, the specific heat transition at high fields shows 
a taller peak in the specific heat than the transition at 
low fields. 

In the magnetocaloric effect data shown in Fig. |2 heat- 
ing is observed as the magnetic field is swept through the 
AFM transition, surrounded by regions of cooling before 
and after the transition. The region of cooling after the 
transition can be attributed to a relaxation towards the 
bath temperature of the system. However, the cooling 
preceding the transition must be due to the sample. The 
inflection point of the T(B) curve, corresponding to the 
point of maximum heat, was identified as the phase tran- 
sition. The data at high field shown in Fig. Efb) shows 
hysteresis between the up and down sweep of the field. 
This hysteresis is observed in the high field branch of 
the phase diagram for temperatures less than 0.6 K, and 
is discussed in greater detail later. The transition tem- 
peratures determined from specific heat and MCE are in 
excellent agreement as shown in Fig. [3] A second order 
phase transition such as the AFM transition we are ex- 
pecting in this compound, should exhibit heating when 
entering the ordered phase, and cooling when leaving it. 
However, all of the MCE data shows heating in both di- 
rections when the field is swept up and down. This may 
indicate a coupling to the lattice that results in an ap- 
parent first-order phase transition. 

The resulting phase diagram (Fig. |3J) shows an ordered 
phase occurring between H c i ~ 2.1 T and H C 2 ~ 12.7 T 
with a maximum critical temperature of T c = 1.2 K. 
The values of H c \ and H c i extracted from this phase 
diagram are in agreement with those determined from 
magnetization to within a few percent^ The discrepancy 
could be due to slight orientation errors of the sample 
and the fact that in the work by Paduan-Filho et al, the 
onset of the transition is determined from the peak in 
the first derivative of M(H), whereas a comparison with 
the midpoint of the transition in the specific heat and 
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FIG. 3: Temperature T - Magnetic field H phase diagram 
from specific heat and magnetocaloric effect data. Bose- 
Einstein Condensation (BEC)/canted antiferromagnetic or- 
der (AFM) is thought to occur in the area under the symbols. 
Top right: Photograph of the single crystal of DTN used in 
this work. The plane of the paper is the a — b plane of the 
crystal. 



MCE effect would require taking the inflection point, or 
the peak in the second derivative of M(H). 

The values of H c \ and for H | | c can be obtained 
as a function of D and J„ using a generalized spin-wave 



approacf 



g^ B H cl = ^D 2 - 2zDJ 



9VbH C 2 = D + 2z J, 



(2) 



(3) 



where the coordination number z = 6. Here J is the 
average of J„ over the three different directions J = 
\{2J a + J c ). The gyromagnetic ratio g c has been deter- 
mined from fits to the magnetization versus temperature 
as g c = 2.26 for H\\c. Thus D and J can be solved for 
yielding D/k B = 10.23 K and J/k B = 0.77 K. In con- 
trast to BaCuSi20e, the DTN compound is in a regime 
for which a two level description does not work because 
2zJ/D <~ 1. In other words, the fluctuations to the high 
energy S z = — 1 state play a role even at very low ener- 
gies. Consequently, the phase diagram is not symmetric 
around the maximum critical field H max (see Fig. [SJl. 
The asymmetry is also manifested as a difference between 
the size of the specific heat jumps at low and high fields 
for a given critical temperature, as shown in Fig. ^ 

In addition to the asymmetry in the shape of the H C (T) 
curve, the MCE data shows hysteresis for T < 0.6 K in 
the high field branch of the phase diagram. Furthermore, 
an anomalous feature can be observed in magnetization 
M(H) data taken by Paduan-Filho el a& at T = 16 mK 
and to a lesser degree at T = 0.6 K. Between H c \ ~ 2 T 
and H ~ 10 T, the M(H) data is linear as expected for 
a BEC. However, above H = 10 T, the slope of M versus 
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FIG. 4: Estimation of the power-law exponent a by fits over 
varying temperature ranges t w = AT/1.2 K (see text), a) Fits 
of the expression H C (T) — H c i oc T a were performed for vari- 
ous fixed values of a at various temperature ranges, resulting 
in Hd a function of temperature range t w . The extrapolated 
value of Hd at zero temperature is H c i — 2.110 ± .002 T. 
b) Using the value of H c i determined in the left figure, fits 
to H C (T) — H c i oc T a axe performed for various temperature 
range to determine a, and a is shown versus t w . The extrap- 
olated value of a at t w = is indicated by the arrow, and the 
lines indicate the error bar. 



H increases until M saturates above H c i- Thus, there 
is a region of anomalous behavior from thermodynamic 
and magnetization measurements between H = 10 T and 
H C 2, and below T = 0.6 K. One possible explanation is a 
gradual expansion in the lattice beginning near H = 10 
T. At high magnetic fields approaching H C 2, the uniform 
magnetization along the c-axis becomes large in compe- 
tition with the antiferromagnetic exchange interaction J. 
Thus, the system can lower its energy by increasing the 
Ni-Ni distance to lower J. As a result, M(H) would 
saturate more rapidly with applied magnetic field. The 
transition from the AFM phase back to the paramagnetic 
phase at high fields would be accompanied by an abrupt 
cessation of the lattice strain, consistent with the hys- 
teretic first-order transition observed in the MCE effect 
at high fields. Further studies, such as magnetostriction 
are necessary to investigate this possibility. 

In order to study the occurrence of BEC of the Ni 
spins in DTN, we examine the temperature dependence 
of the AFM phase transition line near H c \. As men- 
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tioncd previously, a power-law temperature dependence 
H C (T) — H cl oc T a where a = 1.5 is predicted for a transi- 
tion to a BEC phase. In our fits to the H C {T) data for the 
compound DTN, we limit ourselves to investigating the 
low field side of the phase diagram due to the possibility 
of a magnetically induced strain phase occurring at high 
fields. The value of a is closely dependent on the value 
of H c i, as well as on the temperature range of the fit. a 
is expected to approach 1.5 only as T — > 0. The lack of 
particle-hole symmetry does not allow us to use the the 
relation t ~ (1 - h?) v g{h 2 ) {v = 1/a, t = T c /T max and 
h = H max SH 1 ) that was exploited to determine the ex- 
ponent a for BaCuSi2 Therefore, we use the simplest 
power law expression: H C (T) — H c i oc T a . 

As a first step, we fix a and fit the H C {T) data to 
determine H c \. This fit is performed for data up to 
a maximum temperature where the window size t w = 
T max l\.1K . The fits are then repeated for different trial 
values of a. The results of these fits are shown in Fig. 0] 
(left). For each value of a, the H c \ values as a function 
of the window size t w can be fit to a straight line and ex- 
trapolated to a window size t w = 0. The extrapolations 
of H c i(t w ) for the different values of a converge near 
t w = 0, resulting in a mean value of H c i = 2.110 ± 0.002 
T. 

In the expression H C (T) — H c i oc T a , H c \ can now be 
fixed at 2.11 T and the H C (T) data is fit to determine a as 
a function of window size as shown in Fig. 0] (right). For 



t w = 0, a extrapolates to a = 1.47 ±0.06 where the error 
bar in a results from the error bar in H c ±. This value 
is within experimental error of the predicted exponent 
a = 1.5 for a Bose-Einstein condensate. Thus our data 
is consistent with a Bose-Einstein condensation of the Ni 
spin degrees of freedom. 

In summary, we have mapped the high-field phase di- 
agram of DTN using specific heat and magnetocaloric 
effect data, yielding a magnetic ordered state between 
H c i = 2.1 T and H C 2 = 12.7 T with a maximum critical 
temperature of 1.2 K. An anomalous region of the phase 
diagram is observed between H ~ 10 and 12.7 T and 
below T ~ 0.6 K possibly associated with lattice strain. 
The phase transition line near the quantum critical point 
at H c i can be fit by the expression H C (T) — H c i oc T a , 
resulting in a value of a — 1.47±0.06, which is consistent 
with the predicted exponent equal to 1.5 for a quantum 
phase transition to a BEC. 
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